Isospin symmetry breaking at high spins in the mirror pair 67 Se and 67 As 



K. Kaneko 1 , S. Tazaki 2 , T. Mizusaki 3 , Y. Sun 4 ' 5 , M. Hasegawa 4 - 5 , G. de Angelis 6 

1 Department of Physics, Kyushu Sangyo University, Fukuoka 813-8503, Japan 
2 Department of Applied Physics, Fukuoka University, Fukuoka 814-0180, Japan 
^Institute of Natural Sciences, Senshu University, Tokyo 101-8425, Japan 
^Department of Physics, Shanghai Jiao Tong University, Shanghai 200240, People 's Republic of China 
5 Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, People's Republic of China 
^Laboratori Nazionali di Legnaro dell' IN FN, Legnaro (Padova), 1-35020, Italy 
(Dated: November 30, 2010) 

Recent experimental data have revealed large mirror energy differences (MED) between high-spin states 
in the mirror nuclei 67 Se and 67 As, the heaviest pair where MED have been determined so far. The MED 
are generally attributed to the isospin symmetry breaking caused by the Coulomb force and by the isospin 
nonconserving part of the nucleon-nucleon residual interaction. The different contributions of the various terms 
have been extensively studied in the fp shell. By employing large-scale shell model calculations, we show 
that the inclusion of the gg / 2 orbit causes interference between the electromagnetic spin-orbit and the Coulomb 
monopole radial terms at high spin. The large MED are attributed to the aligned proton pair excitations from 
the P3/2 and / 5 y 2 orbits to the gg/ 2 orbit. The relation of the MED to deformation is discussed. 

PACS numbers: 21.10.Sf, 21.30.Fe, 21.60.Cs, 27.50.+e 



One of the current topics in nuclear structure physics is the 
isospin symmetry breaking due to the Coulomb force and the 
strong nucleon-nucleon (NN) interaction. Assuming isospin 
symmetry, mirror pair nuclei, i.e. a pair of nuclei with ex- 
changed proton and neutron numbers, have identical level 
schemes. However, the Coulomb effects and the isospin non- 
conserving NN interaction break this symmetry, leading to ob- 
servable differences between energy levels of analogue states. 
The so-called mirror energy differences (MED) are defined by 

MED/ = E X (J, T, T z = -T) - E X (J, T,T Z = T), (1) 

where E X (J,T,T Z ) are the excitation energies of analogue 
states with spin J and isospin T,T Z . The MED are thus re- 
garded as a measure of isospin symmetry breaking in an effec- 
tive interaction which includes the Coulomb force. The MED 
have been extensively studied for mirror pair nuclei in the up- 
per sd and the lower fp shell regions (see Ref. [ 1] for review). 
In both cases, a remarkable agreement between experimental 
data and shell model calculations has been achieved, allow- 
ing a clear identification of the origin of the MED based on 
the isospin nonconserving Coulomb and strong NN forces Jj- 

a. 

For mirror nuclei in the upper part of the fp shell the situa- 
tion is different. The experimental information on MED is rel- 
atively scarce and only recent large-scale shell model calcula- 
tions including the g 9 j 2 orbit have become available 
Moreover due to the deformation driving effect of the ggj 2 
orbit, variations in the MED are expected to be strongly re- 
lated to the change in the nuclear deformation. Recently new 
data on the A=67 mirror nuclei 67 Se and 67 As have become 
available 112 lh . Investigations B22I1 for 66 Ge suggested that the 
spin alignment of the ggj 2 neutrons occurs at J n = 8+. As 
the positive-parity band built on the 9/2+ state in 67 As can 
be interpreted as a gg/ 2 proton weakly coupled to the 66 Ge 
core, the neutron spin alignment is expected to occur at spin 



25/2+ (= 8+4-9/2+) in 6 'As [22, 23]. On the other hand, the 



proton spin alignment takes place at the same spin in its mirror 
partner 67 Se. As the response to the Coulomb field is different 
for the corresponding high-spin states in such mirror nuclei, 
one expects the Coulomb based MED contribution in 67 Se and 
67 As to give large negative value suddenly at 25/2+ where the 
proton/neutron spin alignment occurs. In the lower fp-shell 
region, due to the active role played by the f 7 1 2 shell, in the 
MED the isospin nonconserving NN interaction has been sug- 
gested to be at least as important as the Coulomb part 12411 . 
In the upper //j-shell region the situation is different and one 
does not expect a major contribution because the f 7 / 2 shell is 
almost not active. 

For the A =67 mirror pair nuclei, excited states have been 
known for 67 As |23|,|25|], and have been recently determined 
for the mirror partner 67 Se ||2g] . This is the heaviest mir- 
ror pair where the excited energy levels have been identified 
with detailed experimental information. In both cases, the 
low-lying 9/2+ state has been found to be isomeric, allowing 
the determination of the degree of isospin symmetry break- 
ing through the measurement of the mirror 9/2+ — > 7 /2~ El 
strengths 112 111 . In our previous paper 119], the structure of this 
isomeric state has been investigated using large-scale shell 
model calculations. The isomerism of the 9/2+ state was 
understood as due to proton and neutron configuration mix- 
ing based on the gg j 2 intruder orbit as well as on the /p-shell 
structures. 

In this Rapid Communication, we investigate the MED in 
the mirror pair 67 Se and 67 As discussing the origin of isospin 
symmetry breaking in the upper fp-shell region. Theoreti- 
cal calculations are performed using the spherical shell model 
in the pf 5/2^9/2 model space. We employ the recently pro- 
posed JUN45 interaction |2(ill . a realistic effective interaction 
based on the Bonn-C potential and ajusted to the experimen- 
tal data of nuclei in the A = 63 ~ 96 mass region. To de- 
scribe the MED, the first attempt was carried out by adding the 
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Coulomb term to the KB3 interaction matrix elements B27I1 . 
However, those calculations did not succeed to describe the 
experimental MED for the mirror pairs of mass A = 47 and 
A =49. A better agreement with the data has been obtained 
using the formalism introduced by Zuker et al. B2411 . In this 
description the Coulomb matrix elements in the valence space 
represent only the multipole part of the Coulomb interaction 
whereas the contribution of the other nucleons is described by 
the Coulomb monopole effect. The Coulomb interaction is 
therefore separated into a monopole term Vc m an d a multipole 
term Vcm- While Vc m accounts for single-particle and bulk 
effects, Vcm contains all the rest. The monopole term Vc m is 
further divided into the single particle correction £//, the radial 
term V& an d the spin orbit term £/ s . The contribution of £// to 
the monopole term is given by f6[] 



-4.5Z 1 J /l2 [2l(l + l)-p(p + 3)] 



AV3(p + 3/2) 



(2) 



where Z cs is the proton number corresponding to a closed 
shell, p the principal quantum number, and / the orbital mo- 
mentum. Due to such single particle correction, in 67 Se the 
proton g 9/ / 2 and / 5 / 2 orbits are lowered roughly by 95 keV 
and 58 keV, respectively, while the energy of the p^/ 2 orbit 
is raised by about 135 keV. The relative energy gap between 
the proton g 9 i 2 and / 5 / 2 orbits is reduced of only 37 keV, and 
therefore there is basically no effect on single-particle levels 
due to the en term. 

The radial term V& reflects the change in radii along the ro- 
tational band, and in the fp shell is proportional to the change 
in occupancy of the p^n orbit as a function of spin /. It can be 
expressed as A M ED{V Cr ) = a m ( (0^3/2)9/2/2 - (m p3 / 2 )j/2), 
where (rtips^j with m p ^/ 2 = Zp3/2 + n p3/2 is the expectation 
value of the proton and neutron number in the 773/2 orbit at 
spin J and a m is the strength parameter fitted to the experimen- 
tal data. When the occupation of the /? 3 / 2 protons decreases, 
valence protons in orbits with smaller radii are nearer to the 
charged core, which results in a gain of Coulomb energy (]]]. 
In the pfe /2g9 /2 shell, the p^ j 2 orbit has larger radius than the 
/s ji and 12 orbits and therefore the Coulomb repulsion in- 
creases as the number of protons increases. Here the role of 
the pi/2 orbit is less important simply because the p X j2 occu- 
pancy is small, and furthermore it does not change very much 
as a function of the angular momentum /. 

The single-particle shift £/, takes into account the relativis- 
tic spin-orbit interaction [28]. This interaction comes from 
the Larmor precession of the nucleons in the electric field due 
to their magnetic moments, which, as well known, affects the 
single-particle energy spectrum, can be written as [28] 
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where m/y is the nucleon mass, and the free values of the gyro- 
magnetic factors, gf=5.586, gf=l for protons and g^=-3.828, 
gj=0 for neutrons, are used. In the present work, by assum- 
ing a uniformly charged sphere, ei s is calculated using the 
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FIG. 1: Comparison of calculated energy levels (JUN45) with exper- 
imental data (exp) for 67 Se and 67 As. 



harmonic oscillator single-particle wave function. Depend- 
ing on proton or neutron orbit, the shift can have opposite 
signs. It depends also on the spin-orbit coupling, as for in- 
stance {l-s) =1/2 when j = 1 + s and (l-s) = — (I + l)/2 
when j = I — s. As this term influences differently on neu- 
trons and protons, its effect becomes very important for some 
particular states. In 67 Se, the proton g 9 / 2 orbit is lowered by 
about 66 keV, while the / 5 / 2 orbit is raised by about 66 keV, 
the effect being opposite for 67 As. Also the relative energy 
gap between the proton gg / 2 and f$ / 2 orbits decreases roughly 
by 132 keV, providing a large contribution to the MED. Since 
the spin-orbit contribution leads to a reduction of the energy 
gap between the proton g 9 / 2 and / 5 / 2 orbits, excitations from 
those orbits into the #9/2 orbit are enhanced. The opposite 
effect is predicted to happen in 67 As for the neutron orbits. 

With inclusion of Vcm, £11 an d £/. s , shell-model calculations 
are carried out in the pf 5 / 2 gt) n shell for the A = 67 mirror 
nuclei. The isospin nonconserving term is neglected in the 
upper half of fp shell region because the / 7 / 2 orbit is almost 
not active. The calculation uses the code MSHELL 12911 and 
the effective interaction JUN45. After solving the eigenvalue 
problem, contribution of the Coulomb monopole radial term 
Vcr is included into the energy Ej obtained in the shell model 
calculation, where the strength parameter a m was fix to 280 
keV so as to fit the MED of the postive-parity high-spin states, 
and taken as 0.0 keV for the negative-parity states. 

In Fig. Q] the calculated energy levels are shown, and com- 
pared with the experimental data for 67 Se and 67 As. As one 
can see, the calculation with the JUN45 interaction reproduces 
well the experimental data. The energy differences of the ana- 
logue states are in a reasonable agreement with experiment. 
The structure of the negative-parity states at low-excitation 
energies are mainly dominated by the fp shell configurations. 
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FIG. 2: (Color online) The MED for states shown in Fig. 1. Upper 
graph: Comparison of calculated MED with available data. Lower 
graph: Decomposition of theoretical MED into four terms (see text 
for explanation). 



The positive-parity states built at higher spin strongly involve 
the orbit. The structural difference of such configura- 
tions strongly reduces the transition strengths explaining the 
isomeric character of the 9/2 + levels IU9I1 . We note however 
that the calculated level energy for the 19/2 + state lies lower 
than the experimental value. 

In the past few years, experimental data on mirror nuclei 
above the doubly magic 56 Ni have become available. Ekman 
et al. Ill 311 discussed the MED of the A ~ 60 mass region based 
on the results of the shell model calculations. It was reported 
that for the MED the contribution of the electromagnetic spin- 
orbit term £/ s is significant, but the monopole Coulomb £// 
term is not. Energy shifts due to Ei s increase the gap between 
the P3/2, /5/2 and the gg/2 orbits for neutrons but reduce it for 
protons. As a consequence excitations involving those orbits 
have important contribution to MED. 

In Fig. |2] (a), the experimental MED along the positive- 
parity excited band with AJ = 2 built on the 9/2 + state and the 
low-lying negative-parity states (3/2~, 5/2~,7/2~) are com- 
pared with the results of our JUN45 calculations as a function 
of spin 27. The agreement is excellent. In particular, the cal- 
culation reproduces correctly the large negative value in the 
MED at the high-spin 21 /2 + and 25/2 + states. It is now in- 
teresting to examine which terms contribute to such drastic 
changes in the MED. In order to see this, the four different 
contributions to MED have been plotted separately in Fig. [2] 
(b). The Coulomb multipole term Vqm reflects the alignment 
effects at high spin and follows the negative trend of the MED. 
It is in fact well known that spin alignments affect the MED, 
which is a behavior first suggested by Sheikh et al. [30] based 
on results from the deformed cranked shell model. In the mir- 
ror pair 49 Mn and 49 Cr the alignment process has been ex- 



amined based on the shell-model calculations by counting the 
number of proton pairs in the shell j coupled to the maximum 
spin J — 2j — 1 [5]. It was shown that this number is closely 
correlated with the MED. For the present case of 67 Se, two 
protons and one neutron jump up from the /p-shell to gg/2 a * 
spin of 25/2+ and 29/2 + . The spin alignment of the two pro- 
tons in the gg / 2 orbit increases the spatial separation between 
them, leading to a smaller Coulomb energy. Thus, the align- 
ment effect for protons reduces the excitation energy in 67 Se 
while the same does not happen in the analogue states in 67 As. 
However, as can be seen in Fig. 12(b), the Vqm term alone un- 
derestimates the MED by a factor of three. As already noticed 
the contribution of the £// term is only marginal. Vo gives the 
largest positive contribution in particular for the 25/2+ and 
29/2 + states due to the increased occupation of the g 9 / 2 orbit. 
On the other hand, the Ei s contribution to the MED is strongly 
negative for the 25/2+ and 29/2+ spin values. When the Ve>, 
Vcm> £/i and en terms are all included, the theoretical MED 
reproduce well the experimental data. We note, however, that 
the strength of Vo was fitted to data and not determined in 
an independent way. The £/ s and V Cr terms contribute to the 
MED from the opposite directions, causing a large cancela- 
tion at the highest spins. For the 21 /2 + state, all terms give 
almost the same contribution of about 40 ke V, providing in to- 
tal large MED. Below 21/2, the Ei s term competes with the £// 
term, while the Vo an d Vqm values are small. In our conclu- 
sion therefore the observed MED behaviour in the 67 Se and 
67 As pair is characterized by a strong competition among the 
different terms, dominated at high spin by the interference of 
the spin orbit and radial contributions. 

A question that we now address concerns the importance 
of the isospin nonconserving term of the NN interaction. The 
good results shown in Fig. [2 have been obtained through the 
inclusion of the Vo term whose strength is however fitted to 
the experimental data. As seen in Fig. |2 calculations without 
Vc r cannot reproduce the data in the high-spin region. In the 
fi /2 shell nuclei the isospin nonconserving NN term is impor- 
tant mainly at the low spin region |l|]. If one speculates that a 
similar behavior occurs also in the g 9 / 2 shell, this would im- 
ply a limited contribution of the isospin nonconserving term 
to the high spin region in the current discussion. The calcu- 
lation presented in this work indeed shows that we can obtain 
a good agreement with the experimental data for the MED 
without including an explicit isospin breaking NN term. All 
these seem to suggest that the isospin nonconserving NN term 
is not important. However, since we have normalised a part of 
the interaction by fitting to the experimental data, we cannot 
make a strong conclusion about the role of the isospin non- 
conserving part that in principle contributes to the MED. 

To support the above picture, Fig. [3] shows the calculated 
occupancies of the excited band with AJ = 2 built on the 9/2 + 
state in 67 Se. The upper and lower graphs are for protons and 
neutrons, respectively. From the upper graph, one can see 
that for the 9/2+ state, protons occupy mainly the /p-shell 
and partially the g 9 / 2 orbit. The occupations change gradu- 
ally such that the fp occupancies increase but the g 9 / 2 one 
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FIG. 3: (Color online) Calculated occupation numbers for proton 
orbits (upper graph) and neutron orbits (lower graph) in 67 Se. 



decreases as a function of 27. However, it is notable that 
the proton gqn occupation increase suddenly at spin 25/2, 
and the proton 773/2 and/5/2 occupations decrease at the same 
spin. This means that two protons and one neutron jump up 
from the //j-shell to the #9/2 orbit at spin 25/2. This drastic 
change of occupations is in clear contrast to that of the / 7 / 2 - 
shell nuclei, where the occupations of P3/2 an d fi/2 orbits 
change gradually with increasing spin yj. The lower graph 
indicates a similar pattern for neutrons, but the variation is not 
so large. As already mentioned above, the change in occu- 
pancy of the /?3 12 orbit affects strongly the MED through the 
Coulomb monopole radial term Vc r - Since the p 3/ / 2 orbit has 
larger radius than the gg / 2 orbit, when at high spin nucleons 
are filling the g<j n shell the Coulomb monopole contribution 
is larger than that at low spins. 

We finally show the calculated spin alignment and spectro- 
scopic quadrupole moment in 67 Se. In Fig. |4] (a), the spin 



distribution of the expectation value J a = y (ja) 2 is plotted 

as a function of spin 27, where j a is angular momentum op- 
erator for each orbit a. As the neutron orbits are blocked for 
this odd-neutron nucleus, the first alignment will be that of a 
pair of g 9 12 protons which brings additional 8 units of angu- 
lar momentum. It is clearly visible that the proton pair and 
one neutron alignment in 67 Se occur at spin 25/2. The 29/2 + 
state also shows a large aligned spin value. This alignment is 
interpreted as five-quasiparticle configuration involving two 
protons and three neutrons. Figure|4](b) shows the calculated 
spectroscopic quadrupole moment Q s (in efm 2 ) for the excited 
states built on the 9/2 + level in 67 Se. The Q s absolute value 
has sudden increase at spin 25/2 corresponding to the sudden 
increase in spin alignment (see the upper graph). This sug- 
gests that the quadrupole moment is closely related to the spin 
alignment of the g 9 / 2 proton pair, which correlates with the 
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FIG. 4: (Color online) Calculated spin distribution in each orbit (up- 
per graph) and spectroscopic quadrupole moment (lower graph) in 
Se. The spin distribution is the spin value J a on a single-particle 
orbit a for the excited states (see text for explanation). 



multipole term Vcm of the MED. Therefore, change in defor- 
mation seems to affect the MED, but its influence is not large. 
It should be noted that the present Vq calculation and the dis- 
cussion on occupation of the single-particle orbits are carried 
out in a spherical basis, and therefore, the deformation effects 
(such as changes of single-particle levels by the shape-driving 
effect) are not explicitly seen. To study the deformation ef- 
fects in the MED, a shell model based on deformed single- 
particle states 1 31 , 32j] would have to be employed. 

In conclusion, we investigated the MED between high-spin 
states in the mirror pair 67 Se and 67 As using large-scale shell 
model calculations. The calculations reproduce well the ex- 
perimental level schemes, and confirm the suitableness of 
the JUN45 effective interaction for this mass region. The 
need for inclusion of the g 9 1 2 orbit in the description for the 
MED in the upper fp shell nuclei was demonstrated. In this 
mass region, the electromagnetic spin-orbit interaction and the 
Coulomb monopole radial term are responsible for producing 
the large MED at high-spin states, while the contribution from 
the Coulomb multipole term is small. The occupations of the 
relevant orbits and the spin alignment in the gg n orbit affect 
the variation of the MED along the band built on the 9/2 + 
state. We obtained a good agreement with the experimental 
data for the MED without involving the isospin nonconserv- 
ing part. However, it cannot be concluded that the isospin 
nonconserving NN term is not important. This remains an 
open question. 
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